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ABSTRACT 


The program is a modification of the Lewis Chemical Equilibrium 
program to extend the usefulness of the program in a Coal 
Gasification Study. A previous version adapted the original 
program to allow direct input of coal composition data. The 
current revision adds a capability for parametric studies of 
various coal gasification factors. In addition, an optional data 
file collection system was implemented for use with the 
parametric studies. The program was also converted to operate 
under ASCII FORTRAN rather than FORTRAN V on the UNIVAC 1100/82. 

The program works on single precision data in an interactive mode 
and uses approximately 32k of core. There are no special 
input/output requirements beyond a single system tape. 
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1.0 INTRODUCTION 
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1.1 OBJECTIVE 

The Coal Gasification Task Team required parametric studies of 
several factors involved in the coal gasification process. The 
parametric studies involve the combination of parameters such as 
oxygen/coal and water/coal ratios over ranges of the individual 
parameters. 

This program was developed to automate the process of running all 
of the possible combinations of the factors under study. The 
factors to be varied are entered as a list of values to be 
computed and the program automatically sets all of the necessary 
combinations of variables and runs the complete gasification 
calculation for each combination. 

The problem which is solved bv the program is the addition of a 
parametric case study capability to the existing Lewis Chemical 
Equilibrium Program. Various parametric data are entered by the 
user and printed data are created covering the possible 
combinations of the original data. 

This program was developed for the Computer Services Office 
(AH33) in support of activities in the Coal Gasification Task 
Team (PF15) under contract NAS8-31540. 

1.2 NSFC FORM 3559 

(See the following page.) 
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1.3 BACKGROUND 


The original Lewis Chemical Equilibrium Program (Reference 1) waa 
developed to determine chemical equilibrium in complex systems. 
Using a free-energy minimization technique, the program permits 
calculations such as: 

e Chemical Equilibrium for assigned thermodynamic states 
(T,P), (H,P) , (S,P), (T,V) , (U,V) , OR (S,V), 

• Theoretical rocket performance for both equilibrium and 
frozen compositions during expansion, 

# Incident and reflected shock properties, 

e Chapman-Jouget detonation properties. The program considers 
condensed species as well as gaseous species. 

tt has been shown (Reference 2) that the same program can handle 
solid coal in an entrained flow coal gasification problem. Since 
the original program was not designed for coal input, the 
calculations necessary to change the coal data into the proper 
form are long and tedious. The automation of the calculation 
process is documented in Reference 2. 

The problem of parametric studies in the original program is 
complicated, because the original program made no provision to 
vary the factors of interest in coal gasification studies. The 
current program was created to reduce the need for user 
intervention at every step of the parametric studies. 


1-3 


2.0 PROGRAM TASK DESCRIPTION 


The Lewis Chemical Equilibrium Program was revised to Include the 
following capabilities* 

1. Accept a data set vftich can be used to create a series of 

runs covering the variations of one parameter as other 

independent variables are changed. All of the variables 
which are being changed affect the coal composition. 
Previous versions of the program had to restate the coal 
composition between each calculation. The current process 
does all of the calculations automatically once the 

calculation values have been entered. 

2. Create, as a user option, a data file which contains a 

condensation of the data calculated in th* multiple runs. 
The option can also be used to collect single data points or 
any combination. 

3. Run under ASCII FORTRAN rather than FORTRAN V. 

A complete listing of the source code, after alterations, is 
given in Appendix A. A complete listing of the original, 

pre- a Iteration, coding is included as Appendix B, so that the 
differences can be traced. 
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3.0 METHOD Of SOLUTION 



I 


The solution of the problem, as described in the Task 
Description, involved creating storage area to hold the 
parametric data and the plot output data. In addition, logic had 
to be added to force the system to loop through the parametric 
values in a planned fashion and to store the plot data as 
required. The conversion to ASCII FORTRAN from FORTRAN V also 
introduced revisions, principally in the area of character 
manipulation and core overlay. 

A detailed description of the effect that the task solution had 
on each subroutine is given in the following pages. 

3.1 STORAGE ELEMENT A 

(Map Overlay Source Code) 

The subroutine COALCV was divided into a driver (COLCV) and 3 
subroutines (C0LCV1, C0LCV2 , C0LCV3 ) t the subroutine 0UT1 was 
divided into 4 subroutines (0UT1, 0UT2, 0UT3, 0UT4) . The 

subroutine SAVE was divided into 2 subroutines (SAVE, NEWOP) . 
The overlay structure was revised to reduce the core requirement, 
since the ASCII level generated larger elements. The split in 
COALCV also ties in with the use of multiple runs, see C0LCV3. 

3.2 STORAGE ELEMENT BLOCK 1 
(Miscellanous data - Block Data Form) 

The Hollerith data was revised to match ASCII : equirements. In 
particular, the contents of COMMON area OUTP were extensively 
revised . 

3.3 STORAGE ELEMENT BOOT 
(Bootstrap Control Section) 

No changes 

3.4 STORAGE ELEMENT CHECK 
(Subroutine Check) 

No changes. Note CHECK is not actually used in the program but 
was included because it is on the system tape. 




3.5 STORAGE ELEMENT COALCV 
(Subroutine COALCV) 

The ordinal COALCV was replaced with a driver for COLCV1, COLCV2 
and C0LCV3. The 3 subroutines, together, are equivalent to the 
old COALCV. The split was arranged partly to produce smaller 
segments but, primarily to enable subroutine C0LCV3 to be entered 
independently of the other two. 

3.8 STORAGE ELEMENT C0LCV1 
(Subroutine C0LCV3) 

The subroutine contains the portion of old subroutine COALCV 
which reads in and displays old coal data, as required,. It also 
determines vdiether the existing coal data is to be completely 
replaced or modified. The particular sections to be modified are 
also identified. 

3.7 STORAGE ELEMENT C0LCV2 
(Subroutine COLCV2 ) 

The subroutine contains the portion of old subroutine, COALCV, 
which reads in the required coal data. This may be an entire set 
of new data or selected replacements of existing data. New 
functions are put in to rearcange the oxidizer list, so that 
oxygen is always first. The ”other fuel” section is rearranged 
so that water is first. Also the variable recycle carbon to coal 
ratio is computed. This value is stored as an additional entry 
in COMMON and CPUEL. It also applies to COALCV, C0LCV1 and 
C0LCV3. 

3.8 STORAGE ELEMENT C0LCV3 

(Subroutine C0LCV3) 

The subroutine contains the portion of the old subroutine COALCV 
which does the calculations to convert raw coal composition data 
into reactant data acceptable to the equilibrium calculation 
program. In addition, the subroutine serves the purpose of 
introducing 3 parameters which are part of the possible 
variations of the multiple runs. 

The three parameters are percent carbon conversion, water to coal 
ratio, and recycle carbon to coal ratio. The parameters enter 
the subroutine as arguments in the call list. If the value of 
the argument is negative, then the existing value of the 
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argument, as taken from the coal data file, is used in the 
calculations. If the argument is positive, then the call 
argument is used in place of the orlqinal argument in all of the 
composition calculations. The subroutine C0LCV3 is called 
separately from CCLCV1 and COLCV2 when changes in the coal 
composition data are required by changes to the three parameters. 
(See subroutine MAIN.) 

The other addition to C0LCV3 t- the storage of certain computed 
date into a storage area of t >e COMMON block DPLOT. Storage 

occurs only if a plot flag has been set in subroutine MAIN. 
Table 1 Indicates the meaning of the various storage values. Not 
all of the storage slots are filled by C0LCV3. Other slots are 
filled by various other subroutines as the data becomes 
available. 


( 


i 
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TABLE 1 


PLOT DATA VALUE8 


VARIABLE 

MUHggR MEANING 


1 

2 

3 

4 

5 
5 

7 

8 

9 

10 
11 
12 
13 


14 

15 
15 

17 

18 

19 

20 
21 
22 
23 


Oxygen/Coal Ratio 

Steam/Coal Ratio 

Carbon Conversion Function 

Recycle Carbon/Coal Ratio 

Fuel/Coal Ratio 

Weight Fraction C in Coal 

Weight Fraction H in Coal 

Weight Fraction N in Coal 

Weight Fraction 0 in Coal 

Weight Fraction S in Coal 

Weight Fraction as in Coal 

Coal Temperature °F 

Higher heating value of Coal Btu/lb 

Steam Temperature °F 

Weight Fraction Oxygen in Oxidizer 

Oxidizer Temperature °F 

Pressure PSIA 

Future Expansion 

Mixture Molecular Weight 
Mixture Temperature °F 
Mixture Enthalpy Btu/Lb 


( 


24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 


Mixture Mole Praction of 
Mixture Mole Fraction of 
Mixture Mole Fraction of 
Mixture Mole Praction of 
Mixture Mole Fraction of 
Mixture Mole Fraction of 
Mixture Mole Fraction of 
Mixture Mole Fraction of 
Mixture Mole Fraction of 
Mixture Mole Fraction of 
Mixture Mole Fraction of 
Mixture Mole Fraction of 
Mixture Mole Fraction of 
Mixture Mole Fraction of 
Mixture Mole Fraction of 
Mixture Mole Fraction of 
liqu ids) 


Future Expansion 



C(S) 

IC(S) 

HCL 

HCN 

nh 3 

S 

(solids & 



3.9 STORAGE ELEMENT CPUS 
(Subroutine CPUS) 


Rearrangement of COMMON area SPECES to allow for more useable 
equivalence statements in other subroutines. Increased array ft 
from 15 x 150 to 25 x 150 to allow for an expansion from 15 to 25 
reactants. This last change was missed in the previous revision. 

3.10 STORAGE ELEMENT DBTON 
(Subroutine DETON) 

Changed COMMON area OUPT as indicated in BL0CK1. Also changed 
Hollerith data to fit into the ASCII character format. 

3.11 STORAGE ELEMENT BQLBRM 

(Subroutine EQLBRM) 

Changed COMMON area SPECES as indicated in Section 3.9. 

3.12 STORAGE ELEMENT FROZEN 
(Subroutine FROZEN) 

Changed COMMON area SPECES as indicated is Section 3.9. 

3.13 STORAGE ELEMENT GAMEFF 

(Subroutine GAMEFF) 

No changes. 

3.14 STORAGE ELEMENT GAUSS 
(Subroutine GAUSS) 

No changes 

3.15 STORAGE ELEMENT GETON 
(Partial control stream) 

The group of control statements labeled GETON is used, via ADD, 
to reinitialize the system after having terminated control. 

3.16 STORAGE ELEMENT HCALC 
(Subroutine HCALC) 

Changed common area SPECES as indicated in Section 3.9. 


3.17 STORAGE ELEMENT LTCPHS 
(Subroutine LTCPHS) 


No changes. 

3.18 STORAGE ELEMENT MAIN 
(System Driver) 

Changes in COMMON area SPECES as indicated in Section 3.9. 
Addition of variable CCR to COMMON area CFUEL. Addition of 
COMMON area DSAVE to provide storage which will not be destroyed 
when subroutine THERMP is overlaid. Addition of COMMON atea 
DP LOT to provide temporary storage for plot data. 

To hold the parameters needed for the multiple runs, the name- 
list CDATA was added. Table 2 shows the components of the name- 
list. Note that for each combination of OXCOAL, PCTCNV, RCCOAL, 
up to six variations of WCOAL may be used. In addition to the 60 
possible combinations which can be used to calculate equilibrium 
compositions, a secondary quench composition can be calculated by 
setting QT to TRUE and assigning values to IQTYPE and TA. A 

total of 120 runs can be created with one set of data with no 
further operator intervention. 

In the initialization section, file 11 is opened via a DEFINE 
FILE. The file is used to accumlate the plot data as one record 
per combination of variables. Also, the material in the CDATA 
list is initialized. 

Immediately after the section which reads in the namelist INPT2, 
a new section is added which does exactly the same thing for 

CDATA. The user is given the choice of using the existing data 
and replacing pieces of it, or initializing the entire group and 

reading in entirely new data. The new CDATA material is stored 

on file 10 after the INPT2 data. 

When execution is actually started, option 7, the user is given a 
choice as to whether or not to print the CDATA material. The 
INPT2 data is always printed. At this point the user is asked 
whether or not plot data storage is required. An answer of 
"yes" requests a code name (16 characters) to associate with the 
data group, sets the plot counter to 1 and than asks whether a 
tape alredy exists. If such a tape exists, then the user should 
already have copied it into file 11. The first record on the 
file gives the number of records on the file and the record 
number of the next record to fill. 
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TABLE 


COAL DATA NAMELIST (CDATA) 


OXCOAL (10) 
PCTCNV (10) 
RCCOAL (10) 
WCOAL (-5,10) 


QT 


IQTYPB 


Oxygen to coal ratio 

Percent carbon conversion 

Recycle carbon to coal ratio 

Water to coal ratio: (Bach 

combination of oxygen to coal , 
percent conversion and recycle 
carbon to coal may have up to 6 
water to coal ratios.) 

True or False, indicates whether a 
quench temperature calculation is to 
be performed. 

Integer which speciifes process 

1 KOPPERS-TOTZEK 

2 TEXACO 


TA 


Quench temperature for Texaco 


The program then computes the number of oxygen to coal 
compositions and the number of sets of water to coal 
compositions. If both are zero, the program will run without any 
use of the multiple run option. 

The temporary parameters of percent conversion, recycle carbon to 
coal and water to coal ratio are set. Note that the maximum 
number of water to coal values is calculated for the combination 
of other variables. The data from INPT2 is obtained from file 10 
and the MIX value is set to the oxygen to coal ratio. When 
quench calculations are required during the second pass (IQT*l)j 
then selected values of the INPT2 variable are set depending on 
which process is being calculated. For the first pass (IQT»0) , 
subroutine C0LCV3 is called to work the new parameters into the 
coal composition data. The processing then continues in the 
normal function, except that after subroutine REACT, the special 
omits and the inert carbon insertions are omitted any time the 
REPEAT flag is on. 

Processing is normal until after the call to subroutine THERMP. 
At this point, the REPEAT flag is set on if any of the flags 
indicating multiple runs are on. Then the programming returns to 
run other cases. The order is Quench calculations, additional 
water to coal ratios and then oxygen to coal ratios. 

When all of the multiple run cases for the set of data have been 
completed, normal processing resumes. Before returning to the 
menu for new input, the stored plot data record on file 11 is 
closed, if the plot flag is on. The number of records associated 
with the given data is added to the first record of the 
particular set. The first record on the file is updated to point 
to the next available record. 

3.19 STORAGE ELEMENT MATRIX 
(Subroutine MATRIX) 

Changed COMMON area SPECES as indicated in Section 3.9. 

3.20 STORAGE ELEMENT NBfOF 
(Subroutine NEW OF) 

The subroutine consists of the entry NEWOF of the old subroutine 
SV/E. Note that there was no overlap between the SAVE and NEWOF 
entries. The only other change was the addition of variable CCR 
to COMMON area CFUEL. 
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3.21 STORAGE ELEMENT OUT1 
(Subroutine OUT1) 

The subroutine consists of the entry OUT1 of the old subroutine 
OUT1 . The old subroutine consists of four entries, OUTl, 0UT2, 
OUT 3 and 0UT4 which were all Independent of each other* The 
entires were split into independent subroutines to reduce core 
requirements in the overlay. 

Other changes were the addition of variable CCR to COMMON area 
CFUEL, the addition of COMMON area DP LOT and the rearrangment of 
some Hollerith data to fit the ASCII requirements. Also, some of 
the plot data is stored. Note that each plot variable has 13 

storage slots. 

During the ordinary multiple run cases, usually only one 
condition is used at a time. The other slots are allocated for 
use when the program is run with multiple pressures or 

tern peratures . 

3.22 STORAGE ELEMENT 0UT2 
(Subroutine 0UT2) 

The subroutine consists of the entry 0UT2 of the old subroutine 
OUTl (see Section 3.20). Other changes were the addition of 

variable CCR to COMMON area CFUEL, the addition of COMMON area 
DPLOT and the revision of COMMON area OUPT, as indicated in 

Section 3.2. Various parts of the plot data are also stored in 
the subroutine. 

3.23 STORAGE ELEMENT OUT3 
(Subroutine 0UT3) 

The subroutine consists of the entry OUT3 of the old subroutine 
OUTl (see Section 3.20). The extra changes are the addition of 
variable CCR to the COMMON area CFUEL and the revision of COMMON 
area OUPT as in Section 3.2 and the revision of COMMON area 

SPECES as in Section 3.9. No plot data is saved in this section. 


3.24 STORAGE ELEMENT OUT4 
(Subroutine 0UT4) 

The subroutine consists of the entry OUT 4 of the old subroutine 
OUTl (see Section 3.20) . Changes in the data areas include the 
Addition of variable CCR to COMMON area CFUEL. the revision of 
COMMON area SPECKS as in Section 3.9 and the addition of COMMON 
area DP LOT. Hollerith data indicating the species which are to 
be stored in the plot area were also added. 

In addition to the original calculations, the mole fractions of 
the various designated species are stored in the plot area. The 
mole fractions of any solids and liquids which are not in the 
designated list are summed and stored in the plot area. 

3.25 STORAGE ELEMENT PLOTBOOT 
(partial Control Stream) 

The set of control statements is used, via ADD. to set the plot 
data file for additions to the file. The + in statement 3 is 
replaced with the actual number of tie tape on which the plot 
data is being accumulated. 

3.26 STORAGE ELEMENT PLOTSAVE 
(Partial Control Stream) 

The control statements are used to transfer the accumulated print 
data from the mass storage file to tape. The tape is the same as 
that from which the data originated. The overwrite feature is to 
be used . 

3.27 STORAGE ELEMENT PROK 
(PDP Element PROK) 

Changed COMMON area SPECES as indicated in Section 3.9. 

3.28 STORAGE ELEMENT REACT 
(Subroutine REACT) 

The changes in the data areas included the addition of variable 
CCR to COMMON area CFUEL. and the addition of COMMON area QUENCH 
and DSAVE . QUENCH contains the material needed to compute the 
quench calculations. DSAVE stores material which must be 
preserved against being overlaid. 
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The only addition to the logic of the subroutine involves the 
quench enthalpy calculation. Quench calculations pass through 
subroutine REACT twice. On the first pass, IQT is ft and normal 
calculations are done. If QT is true, IQT is reset to l and the 
next pass will do the quench calculations. 

The type of calculations depends entirely on the process being 
investigated. For the Koppers-Totzek process, the enthalpy for 
each component of the normal calculation is adjusted by 

Hq-H a - 9.75 X 10" 12 T a 4 (MW) 


Where; 

Hq is the quench enthalpy (cal/g) 

H a is the mixture enthalpy calculated for a normal run 
(eal/g) 

Ta is adiabatic flame temperature corresponding to B a (°K) 

MW molecular weight of the component 

For the Texaco process, the calculation is done by running the 
entire computation at constant T a rather than at the constant 
pressure of the original runs. 

Mote that in either case, when IQT*1, the existing reactant data 
is reused rather than starting from the reactant file data. The 
only difference is that the existing data has all of the 
enthalpies, either calculated or given, while the reactant file 
has only the given enthalpies. For both cases, the asterisk 
which indicates that a given enthalpy should be calculated is 
suppressed in the reactant listing. For the Koppers-Totzek case, 
the enthalpy calculation flag is reset to zero to force a 
recalculation for every component. 

3.29 STORAGE ELEMENT RKTOUT 

(Subroutine RKTOUT) 

Changed COMMON area SPECES as discussed in Section 3.9. Changed 
COMMON area OUPT as discussed in Section 3.2. 
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3.30 STORMS ELEMENT ROCKET 
(Subroutine ROCKET) 


No changes 

3.31 STORMS ELEMENT RREAD 
(Subroutine RREAD) 

The Hollerith data used to form headers and titles was reworked 
to conform to ASCII conventions. 

3.32 STORMS ELEMENT SAVE 
(Subroutine SAVE) 

The subroutine is the portion of the old subroutine SAVE which 
corresponds to the entry print SAVE. (See NEWOP, Section 3.19) . 

3.33 STORAGE ELEMENT SEARCH 
(Subroutine SEARCH) 

Changed COMMON area SPECES as discussed in Section 3.9. 

3.34 STORAGE ELEMENT SHCK 
(Subroutine SHCK) 

Changed COMMON area SPECES as discussed in Section 3.9. Changed 
COMMON area OUPT as discussed in Section 3.2. 

3.35 STORAGE ELEMENT TAPESAVE 
(Partial Control Stream) 

Added statement 3 to get tape numbers of saved tape 
automatically. 

3.36 STORAGE ELEMENT THERMP 
(Subroutine THERMP) 

Changed COMMON area OUPT as discussed in Section 3.2. COMMON 
area DPLOT was added. 

After the call to OUT4, the stored plot data from DPLOT is moved 
to file 11 if the plot flag is on and the record counter is 
incremented by the number of records added. 
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3.37 STORAGE ELEMENT VARPMT 
(Subroutine VARPMT) 

Common area OUPT was changed as discussed in Section 3.2. 

3.38 THERMODYNAMIC DATA (FT LB 4) 

No changes. 

In addition to the previously discussed changes, all 
subroutines were reworked to conform to ASCII level 9 
conventions. Principally, this applied to character data 


of the 
FORTRAN 


4.0 PROGRAM DESCRIPTION 


The program ia run In interactive mode on 
machine under the RXGC 8 operating system 
written in ASCII FORTRAN and requires only 
library routines. 

The program operates in approximately 32K and 
system tape. 


the UNI VAC 1100/82 
. The program ia 
the normal FORTRAN 


uses a single input 
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5.0 OPERATING INSTRUCTIONS 


Table 3 shows a typical set of data for a multiple run. The 
total number of runs to compute this data would be 35 x 2. 
There are 35 water/coal combinations and the quench calculation 
doubles the number of runs. 

The problem involvinq gasification of Kentucky 19 coal in a 
Texaco process using an oxidizer stream of almost pure oxygen and 
with a recycle stream (other fuels) as indicated in Table 3. The 
system pressure is 15 psia and the temperatures of the various 
streams are as indicated. The primary parametric factors under 
consideration are the oxygen to coal ratio and the water to coal 
ratio. Carbon conversion and recycle carbon to coal ratio are 
dependent on the oxygen to coal ratio for this particular coal 
gasification process. 

The following pages trace the steps involved in setting up the 
program and show a typical ^et of output data. 

Note that an answer of ‘YES' to the question “PRINT TO PRINT 
PILE" would place most of the output onto a previously defined 
alternate file (31 ) . 

For option 9. entering the coal data, the particular data already 
existed on the tape. For examples of how to enter a completely 
new set of data or to edit existing data, see "Addition to the 
Lewis Chemical Equilibrium Program to Allow Computation from Coal 
Composition Data". (Reference 2.) 

For option 5, entering namelist variable, the namelist CDATA 
shows the setup to compute all of the requirements of Table 3 in 
one pass. 
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TABLE 3 

TEXACO WITH RECYCLE 


COAL, KY99, T»200°F, P*15 p*i a 

carbon conversion - see table 

other fuels - all 8 T-200°p lb/lb coal 


A1 2 0 3 .07466 

Si0 2 , .1983 

FeO, .08678 

water (1), see table 

Recycle carbon. See table 

0xiiy2er T«22l°F 
.975 0 2 
.025 AR 


Slag 

Nash 


.70 

.65 

.63 

.60 

.55 

.50 

OXYGEN CARBN 
COAL CONV . 

RECYCLE 

CARBN/ 

COAL 



WATER/COAL 



.85 
.90 
.95 
1.00 
1.05 
1 . 10 

.91 37 
.9429 
.9622 
.9750 
.98 34 
.9890 

.0606 

.0379 

.0236 

.0146 

.0086 

.0047 

6087 

5990 

5929 

5890 

5865 

5848 

.7648 

.7526 

.7449 

.7400 

.7368 

.7347 

.8342 

.8209 

.8125 

.8071 

.8037 

.8014 

.9469 

.9318 

.9223 

.9162 

.9123 

.9097 

1.1621 

1.1436 

1.1319 

1.1244 

1.1196 

1.1164 

1. 3977 
1.38 34 
1 . 374 3 
1 . 1684 
1. 364 6 


Ounnch Calculations Required 


Quench Temperature 2000°F 


For option 7, startup, the listing of CDATA shows only the one 
test case which matches the output data. 

Note that an answer of 'YES' to the 'SAVE PLOT DATA' will produce 
a request for a code name of 16 characters to associate with the 
data. It will also produce a question as to whether the plot 
tape is blank or already contains other plot data. 

The first set of output is the normal calculation. The second 
set is for the quench calculation. Note that the special omit 
data section only occurs for the first data set. 

Note that on exiting, if the alternate print file has been 
designated, then the user is required to take the appropriate 
steps to cause the alternate print *ile to be printed. 

5.1 TEST CASE 
See following pages. 
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APPENDIX A 


SOURCE LISTING OP PROGRAM ELEMENTS WITH MODIFICATIONS 

This is a printout of the Lewis Chemical Equilibrium Program with 
the modifications required for the multiple runs. 
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APPENDIX B 


SOURCE LISTING OF PROGRAM ELEMENTS WITHOUT MODIFICATION 


This is printout of the Lewis Chemical Equilibrium Program 
without the modification required for multiple runs. Note that, 
this material is identical to the after modification listing 
(Appendix B) of "Addition to the Lewis Chemical Equilibrium 
Program to «llow Computation from Coal Composition Data." 
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2. Table of Contents 

3. Abstract 

4 Introduction 
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a. Objective or Purpose 

b. MSFC Form 3559 

c. Background 

d. Related 

5. Problem Task Description 

6. Method of Solution 

/. Program Description 

a. Operating Environment 


(1) Hardware 

(a) 

Computer 

(b) 

Core Requirement 

(c) 

Magnetic Tapes 

(d) 

Card Punch 

(e) 

Plotter 
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Drum/Disc 
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Other 

(2) Software 
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Operating System 
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Type of Run 
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Program 
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(4) Magnetic Tape Save Labels . 

(5) Computer Time Requirements 

b. Input 

(1) Cards 

(2) Magnetic Tapes 

(3) Drum/Disc 

(4) Other 

c. Output 

(1) Cards 

(2) Magnetic Tapes 

(3) Drum/Disc 

(4) Printout 

(5) Plots 

(6) Other 

d. Restrictions and/or Limitations 

e. Diagnostics 

f . Test Case 

9. Symbols 

10. References 

11. Appendices 

a. Documentation Checklist 

b. Documentation Approval 
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